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ABSTRACT 



C . We show with several examples that renormalization group (RG) theory can 



be used to understand singular and reductive perturbation methods in a unified 



^ . fashion. Amplitude equations describing slow motion dynamics in nonequilibrium 

H : 

phenomena are RG equations. The renormalized perturbation approach may be 
simpler to use than other approaches, because it does not require the use of 
asymptotic matching, and yields practically superior approximations. 
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The essence of the renormahzation group (RG) is to extract those structurally 
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stable features of a system which are insensitive to details. ' Thus, RG methods 
may be regarded as a means of asymptotic analysis. The usefulness of this point 
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of view has been amply demonstrated by the relation between the RG and in- 

4 

termediate asymptotics, which showed that the anomalous exponents appearing 
in (e.g.) the long-time behavior of certain hydrodynamic systems were calculable 
using RG. 

Many different techniques for asymptotic analysis have been developed in- 
cluding the multiple scaling (MS) method (which actually subsumes all the oth- 
ers), the boundary layer (BL) method, the WKB approximation, and others. 
Reductive perturbation methods have been used to extract the dynamics de- 
scribing the global space-time behavior of complicated systems near bifurcation 
points. 

At a purely technical level, the starting point for both perturbative RG meth- 
ods and conventional asymptotic methods is the removal of divergences from a 
perturbation series. Given the above similarities, a natural question arises: what 
is the relation, if any, between conventional asymptotic methods and the RG? 

In this Letter, we demonstrate that many singular perturbation methods 
may be understood as renormalized perturbation theory, and that amplitude 
equations obtainable by the reductive perturbation methods are renormahzation 

o 

group equations. One of the advantages of the RG approach is that the starting 
point is simply a straightforward naive perturbation expansion, for which very 
little a priori knowledge is required; we will see that the RG approach seems to 
be more efficient and accurate, in practice, than standard methods in extracting 
global information from the perturbation expansion. 



To illustrate the basic idea, let us consider a weakly nonlinear van der Pol 
oscillator 

d^y , (dy 1 fdyYi , . 

5 
which is usually solved by MS. The method of uniformization or 'renormalization' 

fails in this case. A naive expansion y = yo + eyi + <P'y2 + ■ ■ ■ gives 

?2~ 



y{t) =Ro sin(t + Go) + e{ Aii?o sin(t + Gq) + (bi - ^ j Rq cos(t + Gc 
+ ^(l-^Vt-to)sin(t + eo) + ^cos3(t + eo)}+0[e2] 



(2) 



where Rq^Qq, Ai,Bi are constants determined by the initial conditions at arbi- 
trary time t = tQ. This naive perturbation theory breaks down for e(t — to) > 1 
because of the secular terms. The arbitrary time to J^ay be interpreted as 
the (logarithm of the) ultraviolet cutoff in the usual field theory. To regular- 
ize the perturbation series, we introduce an arbitrary time r, split t — to as 
t — r + r — to, and absorb the terms containing r — to into the renormalized 
counterparts R and of Rq and Go, respectively. This is allowed since Rq and 
Oo are no longer constants of motion in the presence of the nonlinear perturba- 
tion. We introduce a multiplicative renormalization constant Zi = 1 + YlT ^n^^ 
and an additive one Z2 = YlT ^n^^ such that -Ro(^o) = -^1(^07 ''")-R(t) and 

©oli^o) = ^{t) + ■^2(to,T). The coefficients an and bn {n > 1) are chosen or- 

9 
der by order in e to eliminate the terms containing r — to as in the standard RG. 



The choice ai = — (1/2)(1 — i?^/4)(r — to) — Ai, 61 = —Si removes the secular 
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terms to order e, and we obtain the following renormalized perturbation result 

y{t) =\r + e{R/2){l - R^/A) (t - r) | sin(t + 6) 

^ J (3) 

- e(l/96)i?3cos(t + e) + e(i?V96) cos3(t + 6) + 0[e\ 



where R, G are now functions of r. Since r does not appear in the original 
problem, the solution should not depend on r. Therefore, {dy/dT)t = 0. This is 
the RG equation, which in this case consists of two independent equations 

dR/dr = e{R/2) (l - R^/4) + 0[e^], dG/dt = 0[e^]. (4) 

Setting r = t in (3) eliminates the secular term, giving 

y{t) = R{t) sin(t) + {e/96)R{tf cos(3t) + 0[e^], (5) 

where R{t) and 6(t) are obtained from (4) with t = t: R{t) = R{0)[e^'^^ + 
i?(0)2(l-e~^*)/4]"^/2 and e(t) = 0(0), which we take to be zero. The final result 
approaches a limit circle of radius tending towards 2 as t ^ oo. The second order 
RG calculation shows that our assumption of perturbative renormalizability is 
consistent and no ambiguity arises. The multiple time scales used in MS for slow 
variables Ti = et, T2 = e^t, ■ ■ ■ appear naturally. 

The above example illustrates two important points: (1) the results of the MS 
method can be obtained from renormalized perturbation theory, and (2) the RG 
equation describes the long time scale motion of the amplitude and the phase. 
In the following we wish to demonstrate (1) and (2) more generally. 

Another important class of singular problems is that for which the highest 

order derivative of the equation is multiplied by a small parameter e, e.g. WKB 

5 
and BL problems. For linear cases, it is known that both problems can be 

treated in a unified fashion. A typical problem is of the form 

'—^ + a{x)-^-b{x)y = 0, < x < 1, e ^ 0+, (6) 
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where a and b are continuous functions, and we have chosen a{x) > so that 



the boundary layer is at x = 0. This can be transformed into a form suitable for 
WKB analysis: 



dx"^ 



Q{x)u{x), 



(7) 



with Q{x) = a^(a;)/4e^ + a'{x)/2 + b{x) and 



y{x) = exp 



X 

-—^ / C'i^')dx' 



u{x) 



(8) 



It is convenient to introduce a new variable t such that dt — ^fQdxje^ and (7) 
becomes 



dt^ ^-~\Q'"'^^^))Q'^^^^)% 



(9) 



where Q'{x{t)) is assumed to be a function of order unity varying slowly on the 
time scale t, and Q{x) ^ holds for < x < 1 for simplicity. Naively expanding 
u as u{i) = uo{t) + eui{t) + ■ ■ ■, we find the bare perturbation result 



u{t) =e*|Ao + eAo f S {x{s)) ds - eAoe'^^ f S {x{s)) e^'ds\ 



to 



to 



+ c'^IBq + eBo f S {x{s)) ds - eSoe^* f S {x{s)) e-^'ds\ + 0[e^ 



to 



to 



(10) 



where S{x) = —Q^^''^Q'{x)/A, and Aq, Bq are constants dependent on the initial 
conditions at t = to- The terms in the curly brackets are the secular terms 
divergent in the limit t — to -^ oo. These terms are renormalized away with 
the aid of the multiplicative renormalization Aq = ZiA{t) and Bq = Z2-B(r), 
where A, B are the renormalized counterparts of Aq, Bq, respectively. Here Zi — 



1 + ej ° S {x{s)) ds + ■■■ = ^27 with r being some arbitrary time, as in the 
example of the van der Pol oscillator. The renormalized perturbation result is 



u 



(t) = e^{A(r) + eA(r) ( Sds\ + e-^{s(r) + eS(r) ( Sds\ + 0[e], (11) 



where 0[e\ refers to all regular terms of order e which remain finite even as 
t — to ^ oo. The RG equation du/dr = gives 



'^^+e\Q-^I^Q'{x{r))C = 0[e\ 



(12) 



where C = Aox B. Again, (12) corresponds to the amplitude equation with t = t, 
which gives C{x) ~ Q^^'^{x). This is just the adiabatic invariant A{x)Q^'^{x) — 
A(0)(5^/^(0) = constant. Therefore, the 'physical-optics' approximation for (7) 
reads 



u{x) ~ Q-^'\x) { Ci 



exp 



dx' \/ Q{x' 



+ C2 exp 



dx' \IQ{x' 



(13) 



where C\ and Ci are integration constants. The asymptotic result y{x) uniformly 
valid for ex <\ for the general linear boundary-layer problem (6) is given by (8) 
with an appropriate asymptotic expansion formula for Q{x). Notice that the 
above RG approach gives uniformly reliable results from the "inner expansion" 
alone without any (intermediate asymptotic) matching. 

We have seen that the RG equation becomes the equation of motion for the 
slow behavior of the system. To see how RG reproduces reductive perturbation 
results, let us consider the following equation (this type covers most examples so 



far studied in the literature) 

[Ci{dt)+C2{V)]u = eQ[u], (14) 

where £i and £2 are constant coefficient hnear differential operators, Q is a 
possibly nonlinear operator, and e is a small parameter. We assume, for simplic- 
ity, spatial isotropy. Suppose that the operators have the following structures 

^i = Yluidt + zu;)"^("), £2 = n^(V - iy)<^\ and uq = E.;,^ ac^.^e^^^'"""') is a 
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solution to (14) with e = 0. The order e correction, -ui, in the naive perturbation 

obeys (£1 + C2)ui = Q[uq]. We assume without any loss of generality that Q[uq] 

can be expanded as 

Q[uo] = Yl Qc.,;.[{a}]e^('^-^-^^) + R, (15) 

LJ, fj, 

where R is the remainder such that CiR 7^ 0, C2R 7^ 0, and Q^;,^ are coefficients 
dependent on the set of atj,^'s collectively denoted by {a}. The general form of 
the singular (secular) part [ui]s of the general solution is 

(16) 

+ (1 - A)£2(^)-^ '^' ^^^^7' + p.,, e^(---*) 



; 



where A is an arbitrary numerical constant, not equal to or 1, P^,^ is a 
polynomial of t and |xp of lower order than m{uj) and n(/u)/2, respectively 
(whose explicit form is not required), i-i{uj) = Ylcj'^uji''''^ ~ iuj')'^^'^ ' and i-2{p') — 
Yla'^ui'^f^' ~'^f^) ■ Renormalization of the secular terms divergent in the global 



space-time limit can be done following the procedures given above, and is tan- 
tamount to replacing in (16) t" — Iq with t" — r", |xp with \x\'^ — \r\'^ and 
the 'bare' coefficients {a} with their renormalized counterparts {A}, regarded as 
functions of r and r. The renormalization group equation can be obtained from 
the condition that u is independent of the parameters r and r introduced in the 
renormalization process. The term P is dependent upon e.g., initial conditions. 
Thus, to obtain a universal result, we differentiate u sufficiently many times with 
respect to r and r to eliminate P. Further eliminating A, we find the following 
mode-coupled amplitude equation: 

^i(^) g^n^H'^^ +^sMA^^^^/^A^,^ = eQ^AlA}). (17) 

Here we have used the isotropy to introduce the Laplacian A. This is the renor- 
malization group equation independent of the arbitrariness due to initial condi- 
tions, solution methods, etc. 
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As an example, consider the Swift-Hohenberg equation: 

du/dt^eu-{l + V'^)^u-u^, (18) 

where e is a control parameter. Here £i = dt, £2 = (1 + V^)^, eQ[u] = eu — u^ 

and uo = ae*^ + a*e-^^. Thus £1 = 1, i2i±l) = -4, eQo,i {A A*) = eA - 3\A\'^A 
and eQo-i{A,A*) = eA* - 3\A\'^A*. That is, 

^ = 4g^ + M-3|-4p-4). (19) 

Scaling out e identifies the slow time and spatial variables asT = et,X = xjl^. 
The scaled result is the well-known time-dependent Ginzburg-Landau amplitude 



equation. Because the secular terms et and ex"^ have been removed from (16) the 
outcome (19) should be uniformly valid up to time scale 1/e and spatial scale 
1/v^ for e ^ 1. In our approach, spatial and time coordinates are treated on an 
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equal footing, in contrast to the standard reductive perturbation method. 

To demonstrate that there are not only conceptual but also technical advan- 
tages to the RG approach, we conclude with a problem involving the so-called 
'switchback': conventionally, only through subtle analysis in the course of actu- 
ally solving the problem is it possible to realize the need for (e.g.) unexpected 
order terms to make asymptotic matching consistent. An example is a carica- 
ture of the Stokes-Oseen singular boundary layer problem, which describes the 
low Reynolds number viscous flow past a sphere of unit radius. After scaling 
the radial coordinate r to a; = er, where e is the Reynolds number squared, the 
equation is 

(Pu 2 du du , , , , , , 

— ^ H — + u—- = 0, u(x = e) = 0, u(x = oo) = 1. (20) 

aa;^ x dx dx 

We regard (20) as an initial-value problem, given an initial condition u{xq) = Aq 
at some arbitrary point x = xq, where Aq is a flnite constant. Assuming a naive 
expansion u{x; e) = uq{x) + Xi{e)ui{x) +A2 (6)^2(2;) + ■ ■ ■, where Aj(e), z = 1, 2, ■ ■ • 
will be determined in a self-consistent way, we obtain uq{x) = Aq. Solving to 
0[Ai(e)] gives 

u{x) ^Aq + Xi{e)AiAo{e2{AoXQ) - e2{Aox)} + 0[Xi{e)^, A2(e)], (21) 

where Ai{xq) is some constant of integration, and we define e2(t) = 
/^°° dpp~^e-P ~ 1/t + logt + (7 - 1) - t/2 + 0{f) ast^O with Euler's con- 
stant 7 ~ 0.577. The naive perturbation result (21) breaks down when both 



xq is small and x — xq is large. To cure this we introduce the renormalization 
constant Z = 1 — X{e)Ai{e2{An) — 62(^x0)} such that A{fx) = ZAq, giving the 
renormalized perturbation series 



u 



(x) = A{ij) + Xi{e)AiA{e2{Aij) - 62(^0;)} + 0[Ai(e)2, A2(e)]. (22) 



The RG equation du/dfx — yields, after setting n — x, 

_Li = Xi{e)Ai^-^ + 0[Ai(e)2, A2(e)], (23) 

and we find the final uniformly valid asymptotic result u{x) = A(x). 

Equation (23) can be solved iteratively along with the required boundary 
conditions A{oo) = 1 and A{e) = as 

CX) 

A{x) = l-Xi{e)Ai fdpp-\-P + 0[Xi{e)^X2{e)]. (24) 

X 

The condition A{e) = gives Ai(e)Aie2(e) = 1. Therefore, the expansion co- 
efficient Ai can be explicitly chosen as Xi{e)Ai = l/e2(e), with the asymptotic 
expansion, in the limit e -^ 0+, AiAi(e) ~ e(l — eloge — (7 — l)e + ■•■). In 
addition, we require that Ai /A2 = 0(1), so that the equation for U2 yields new 
information. The resulting asymptotic solution is correct to 0[e loge] and agrees 
with that obtained by asymptotic matching. Note that in our method, the eloge 
term appears naturally from the asymptotic expansion of 62 (e), whereas some 
artistry is required to obtain this term conventionally. The result to 0[e] given 
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by asymptotic matching is obtained from the renormalized perturbation ex- 
pansion to 0[A2]. The asymptotic expansion to 0[e] is not uniformly valid in r, 
and a much better approximation, in practice, is our full result 1 — e2(er)/e2(e) 
to order Ai. This can be seen clearly in figure 1. 
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In summary, we have demonstrated that various singular perturbation meth- 
ods and reductive perturbation methods may be understood in a unified fashion 
from the renormahzation group point of view, with some attendant technical 
advantages. 
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FIGURE CAPTIONS 

Fig. 1. Comparison between the numerical solution of eq. (20) for several values of 
e, the first order RG result 1 — e2(er)/e2(e), and two matched asymptotic 
expansions (one at fixed r, the other at fixed p = re), as derived in ref. 14. 
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